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Abstract. A new determination of the pregalactic helium abundance 
based on the Magellanic Clouds H II regions is discussed. This deter- 
mination amounts to Y p = 0.2345 ± 0.0030 and is compared with those 
derived from giant extragalactic H II regions in systems with extremely 
low heavy elements content. It is suggested that the higher primordial 
value derived by other authors from giant H II region complexes could 
be due to two systematic effects: the presence of neutral hydrogen inside 
the helium Strdmgren sphere and the presence of temperature variations 
inside the observed volume. 



1. Introduction 

The determination of the pregalactic, or primordial, helium abundance by mass 
Y p is paramount for the study of cosmology, the physics of elementary particles, 
and the chemical evolution of galaxies (e. g. Fields & Olive 1998, Izotov et 
al. 1999, Peimbert & Torres-Peimbert 1999, and references therein). In this 
review we briefly discuss the method used to derive Y p and its main sources of 
error as well as a new determination based on observations of the SMC. This 
determination is compared with those carried out earlier based on extremely 
metal poor extragalactic H II regions. 

The Magellanic Clouds determination of Y p can have at least four significant 
advantages and one disadvantage with respect to those based on distant H II 
region complexes: a) no underlying absorption correction for the helium lines 
is needed because the ionizing stars can be excluded from the observing slit, b) 
the determination of the helium ionization correction factor can be estimated 
by observing different lines of sight of a given H II region, c) the accuracy 
of the determination can be estimated by comparing the results derived from 
different points in a given H II region, d) the electron temperature is generally 
smaller than those of metal poorer H II regions reducing the effect of collisional 
excitation from the metastable 2 3 S level of He I, and e) the disadvantage is that 
the correction due to the chemical evolution of the SMC is in general larger than 
for the other systems. 
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2. He + /H + Determinations 

To derive accurate He + /H + values we need very accurate N e (He II), T e (He II), 
and r(3889, He I) values. Good approximations to determine He + /H + for N e < 
300 cm" 3 , 13 000 K < T e < 20 000 K and r(3889) = 0.0, have been presented 
for the main helium lines by Benjamin, Skillman and Smits (1999): 



iV(He+) 1(6678) n ^0.249-2.0x10-^, 



N(H+) J(H/3) 



2.58T 4 U - Z4S ™ X1U J \ (1) 



7(4471) nl ^Q,i27-4.1xlQ- 4 Af e / \ 

Tim 4 ' [) 

1(5876) „ QC .„o.230-6.3xlO- 4 Af e /q-v 

7^0-735i 4 . (3) 



Fortunately for giant extragalacic H II regions r(3889) is very small and 
frequently close to zero. r(3889) can be estimated together with iV e (He II) 
from the 3889/4471 and the 7065/4471 ratios computed by Robbins (1968). 

Most authors assume that T e (0 III) = T e (He II), there are two reasons for 
this assumption the He + and ++ emission regions occupy similar volumes (but 
not identical ones) and T e (0 III) is easy to measure. Nevertheless the assumption 
is correct only for isothermal nebulae, and not for real nebulae if very accurate 
abundances are needed. In the presence of temperature variations along the line 
of sight the [O III] lines originate preferentially in the high temperature zones 
and the helium and hydrogen lines in the low temperature zones (e.g. Peimbert 
1967, 1995). It can be shown that the temperature derived from the ratio of the 
Balmer continuum to a Balmer emission line, T e (Bac), is similar to T e (He II), 
and that both temperatures are smaller than T e (0 III). 

From models computed with CLOUDY (Ferland 1996) it is found that 
T e (Bac) —labeled T e (Hth) by CLOUDY— is about 5% smaller than T e (0 III). 
Moreover in the model with a homogeneous sphere of I Zw 18 by Stasinska and 
Schaerer (1999) it is found that (T e (Ar III)} = 16300 K (for 63% of the volume) 
and (T e (Ar IV)) = 18300 K (for 36% of the volume) indicating the presence of 
temperature variations. Notice that the average temperature has to be weighted 
by the emissivities strengthening the effect of the temperature variations. More- 
over there is additional evidence that indicates that the temperature variations 
are even higher than those predicted by photoionization models (e. g. Peim- 
bert 1995; Luridiana, Peimbert, & Leitherer 1999; Stasinska and Schaerer 1999). 
From equations (1-3) it follows that the smaller the adopted temperature the 
smaller the derived He + /H + value. 

Due to the presence of very strong density variations inside gaseous nebulae 
the different methods to derive the density yield very different values. The root 
mean square density, iV e (rms), is usually obtained from the observed flux in a 
Balmer line and by assuming a spherical geometry; iV e (rms) provides a minimum 
value for N e (He II), the local density needed to derive the helium abundance. 
Forbidden line ratios of lines of similar excitation energy give us an average 
density for cases where the density is similar or smaller than the critical density 
for collisional deexcitation; available line ratios in the visual region are those of 
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[S II] , [O II] , [CI III] , and [Ar IV] , unfortunately giant extragalactic H II regions 
are close to the low density limit of these ratios, the line intensities of [CI III] 
are very faint, for most observations of the [O II] lines the resolution is not high 
enough to separate A3726 from A3729 nor A4711 of [Ar IV] from A4713 of He I. 
Consequently most of the densities in the literature are those derived from the 
[S II] lines, but as rightly mentioned by Izotov, Thuan, & Lipovetsky (1994, 
1997) they are not representative of the regions where the He I lines originate. 
The self-consistent method, advocated by Izotov and collaborators, is based only 
on line ratios of helium I and is the best method to derive the density of the 
He II zone: N e (Re II). 

The stellar underlying absorption can affect the derived He I and H I emis- 
sion line intensities and has to be taken into account. The best way to reduce 
this effect is to avoid the presence of bright early type stars in the observed 
slit, this can only be done for objects inside the Galaxy and the Magellanic 
Clouds. To minimize the effect of the underlying absorption it is recommended 
to use only objects where the line intensities show very large equivalent widths 
in emission, and to increase the spectral resolution. The best way to correct 
for underlying absorption is to use starburst models that predict the underly- 
ing stellar spectrum. The correction for underlying absorption can be tested by 
comparing the higher order Balmer lines, that are most affected by this effect, 
with the brightest Balmer lines that are the least affected. Similarly the un- 
derlying absorption effect is larger for A4471 and smaller for AA5876, 6678 and 
7065; the effect for these three lines is in general negligible due to a combination 
of causes (mainly their large equivalent widths in emission). 



3. Ionization Structure 

The total He/H value is given by: 



AT (He) AT (He ) + V(He+) + iV(He ++ ) 

A/p)" " AT(H°) + AT(H+) ' ^ ' 

= CT( He,«±«. (5) 

The He ++ /H + ratio can be obtained directly from the 4686/H/3 intensity 
ratio. In objects of low degree of ionization the presence of neutral helium inside 
the H II region is important and ICF(He) becomes larger than 1. The ICF(He 
can be estimated by observing a given nebula at different lines of sight since He' 
is expected to be located in the outer regions. Another way to deal with this 
problem is to observe H II regions of high degree of ionization where the He 
amount is expected to be negligible. Vflchez & Pagel (1988) (see also Pagel et 
al. 1992) defined a radiation softness parameter given by 

_ N(Q + )N(S ++ ) _ 
' N(S+)N(0++) , {) 

for large values of n the amount of neutral helium is significant, while for low 
values of n it is negligible. 
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Figure 1. Spectra of NGC 346 with and without underlying absorp- 
tion. The vertical scale is for the lower spectrum (region 346-A). The 
flux of the upper spectrum (region 346-B) was normalized to the Ha 
emission line flux of the lower spectrum. 

On the other hand for ionization bounded objects of very high degree of 
ionization the amount of H° inside the He + Stromgren sphere becomes signifi- 
cant and the ICF(He) can become smaller than 1. This possibility was firstly 
mentioned by Shields (1974) and studied extensively by Armour et al. (1999) 
for constant density models. Since H° will be located in a thin shell in the 
border of the nebula it can be shown that for ionization bounded nebulae of 
homogeneous density the fraction of H°/He + will be three times larger for an 
observation of the whole object than for an observation of a line of sight that 
includes the center. Similarly this fraction becomes higher than a factor of three 
for models with decreasing density from the center outwards. Alternatively for 
density bounded nebulae this effect can be neglected. 



4. SMC 

Peimbert, Peimbert, & Ruiz (2000a) presented long slit observations of the most 
luminous H II region in the SMC: NGC 346. They divided the two long slit 
positions into thirteen areas, four including the brightest stars {m ~ 14) and 9 
without stars brighter than m = 17. In the upper part of Figure 1 we present a 
spectrum that includes all the observed areas (region 346-B); while in the lower 
part of Figure 1 we present a spectrum made by the seven brightest areas that 
do not include the brightest stars (region 346-A). 

After correcting region 346-A for extinction, based on the four brightest 
Balmer lines, it is found that the weaker Balmer lines (H9 to H12) are not 
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affected by stellar underlying absorption (see Figure 1), therefore the He lines 
are not expected to be affected by underlying absorption. The extinction derived 
from the Balmer lines is in very good agreement with the extinction derived from 
the stellar cluster by Massey, Parker, & Garmany (1989), another indication that 
the stellar underlying absorption of the Balmer lines is negligible. 
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a Given in units of 10 4 , x 2 values in parenthesis. 

b The He + /H + values in boldface correspond to the minimum x 2 values at 
a given temperature. 

c The minimum x 2 value at a given temperature is presented in italics. 

d The smallest x 2 value for all temperatures and densities is presented in 
boldface, thus defining T e (He II) and 7V e (He II). 



To derive the He + /H + value, in addition to the Balmer lines, we made use 
of nine He I lines, AA 3889, 4026, 4387, 4471, 4921, 5876, 6678, 7065, and 7281 
to determine A r e (He II) and T e (He II) self-consistently. In Table 1 we present 
He + /H + values for different temperatures and densities; the temperatures were 
selected to include T e (0 III), T e (Bac), T e (He II), and two representative temper- 
atures; the densities were selected to include the minimum x 2 at each one of the 
five temperatures. The temperature with the minimum x 2 1S the self-consistent 
T e (He II) and amounts to 11950 ± 560 K; this temperature is in excellent agree- 
ment with the temperature derived from the Balmer continuum that amounts 
to 11800 ± 500 K, alternatively T e (0 III) amounts to 13070 ± 100 K. Notice 
that the x 2 test requires a higher density for a lower temperature, increasing 
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the dependence on the temperature of the He + /H + ratio. The values in Table 1 
correspond to the case where r (3889) equals zero, for higher values of r (3889) 
the x 2 values increase. In Table 1 the He + /H + values in boldface and the % 2 
values in italics correspond to the minimum value of x 2 a t a given T e , the x 2 
value in boldface is the smallest value for all temperatures and all densities. 

From Table 1 we obtain that T e (He II)= 11950 K and iV e (He II)= 143 cm" 3 , 
which correspond to He + /H + = 0.0793 ± 0.007. By comparing the He/H values 
for lines of sight with different ionization degree it is found that ICF(He) = 1.00. 
To obtain the total He/H value we have added the contribution of He ++ /H + 
that amounts to 2.2 x 10 -4 . 

In Table 2 we present the helium abundance by mass Y(SMC), derived from 
NGC 346. The Y(SMC) values were derived from Table 1, the He++/H+ values 
and the Z value determined by Peimbert et al. (2000a). 



Table 2. Y(SMC) 

T e (K) \. 3 (cm- 3 ) 

53 100 143 162 247 



11200 
11800 
11950 
12400 
13000 



0.2431 
0.2435 
0.2436 
0.2439 
0.2443 



0.2416 
0.2419 
0.2420 
0.2421 

0.2423 



0.2404 
0.2405 
0.2405 

0.2406 
0.2407 



0.2399 
0.2399 

0.2399 
0.2399 
0.2400 



0.2377 a 

0.2375 
0.2374 
0.2372 
0.2370 



a Boldface values correspond to minimum x 2 values, see Table 1. 



To determine the Y p value from the SMC it is necessary to estimate the frac- 
tion of helium present in the interstellar medium produced by galactic chemical 
evolution. We will assume that 

Y p = Y(SMC) - Z(SMC)^|. (7) 

Peimbert et al. (2000a) find that for T e (He II)= 11950 ± 500 K, Y(SMC)= 
0.2405 ± 0.0018 and Z(SMC)= 0.00315 ± 0.00063. To estimate AY/AZ we will 
consider three observational determinations and a few determinations predicted 
by chemical evolution models. 

Peimbert, Torres-Peimbert, & Ruiz (1992) and Esteban et al. (1999) found 
that Y = 0.2797±0.006 and Z = 0.0212±0.003 for the Galactic H II region M17, 
where we have added O.lOdex and 0.08dex to the carbon and oxygen gaseous 
abundances to take into account the fraction of these elements embedded in 
dust grains (Esteban et al. 1998). By comparing the Y and Z values of M17 
with those of NGC 346 (Peimbert et al. 2000a) we obtain AY/AZ = 2.17 ± 
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0.4. M17 is the best H II region to determine the helium abundance because 
among the brightest galactic H II regions it is the one with the highest degree 
of ionization and consequently with the smallest correction for the presence of 
He (i.e. ICF(Re) is very close to unity). It can be argued that the M17 values 
are not representative of irregular galaxies, on the other hand they provide the 
most accurate observational determination. From a group of 10 irregular and 
blue compact galaxies, that includes the LMC and the SMC, Carigi et al. (1995) 
found AY/ AZ = 2.4±0.6, where they added 0.2dex to the O/H abundance ratios 
derived from the nebular data to take into account the temperature structure 
of the H II regions and the fraction of O embedded in dust, moreover they also 
estimated that O constitutes 54% of the Z value. Izotov & Thuan (1998) from 
a group of 45 supergiant H II regions of low metallicity derived a AY j AZ = 
2.3 ± 1.0; we find from their data that AY j AZ = 1.5 ± 0.6 by adding 0.2dex 
to the O abundances to take into account the temperature structure of the H II 
regions and the fraction of O embedded in dust. 

Based on their two-infall model for the chemical evolution of the Galaxy 
Chiappini, Matteucci, & Gratton (1997) find AY I AZ = 1.6 for the solar vicinity. 
Carigi (2000) computed chemical evolution models for the Galactic disk, under 
an inside-out formation scenario, based on different combinations of seven sets of 
stellar yields by different authors; the AY j AZ spread predicted by her models 
is in the 1.2 to 1.9 range for the Galactocentric distance of M17 (5.9 kpc). 

Carigi et al. (1995), based on yields by Maeder (1992), computed closed box 
models adequate for irregular galaxies, like the SMC, and obtained AY I AZ = 
1.52. They also computed models with galactic outflows of well mixed material, 
that yielded AY j AZ values similar to those of the closed box models, and 
models with galactic outflows of O-rich material that yielded values higher than 
1.52. The maximum AY j AZ value that can be obtained with models of O-rich 
outflows, without entering into contradiction with the C/O and (Z—C — 0)/0 
observational constraints, amounts to 1.69. 

Carigi, Coh'n, & Peimbert, (1999), based on yields by Woosley, Langer, 
& Weaver (1993) and Woosley & Weaver (1995), computed chemical evolution 
models for irregular galaxies also, like the SMC, and found very similar values 
for closed box models with bursting star formation and constant star formation 
rates that amounted to AY I AZ = 1.71. The models with O-rich outflows can 
increase the AY/AZ, but they predict higher C/O ratios than observed. 

From the previous discussion it follows that AY/AZ = 1.9 ± 0.5 is a repre- 
sentative value for models and observations of irregular galaxies. Moreover, this 
value is in good agreement with the models and observed values of the disk of 
the Galaxy. 

The Y p values in Table 3 were computed by adopting AY/AZ = 1.9 ± 0.5 . 
The differences between Tables 2 and 3 depend on T e because the lower the T e 
value the higher the Z value for the SMC. 

6. Discussion 

The Y p value derived by us is significantly smaller than the value derived by 
Izotov &: Thuan (1998) from the Y - O/H linear regression for a sample of 45 
BCGs, and by Izotov et al. (1999) from the average for the two most metal 
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Table 3. Y p Derived from the SMC 
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a Boldface values correspond to minimum \ 2 values, see Table 1. 



deficient galaxies known (I Zw 18 and SBS 0335-052), that amount to 0.2443 ± 
0.0015 and 0.2452 ± 0.0015 respectively. 

The difference could be due to systematic effects in the abundance determi- 
nations. There are two systematic effects not considered by Izotov and collabo- 
rators that we did take into account, the presence of H° inside the He + region 
and the use of a lower temperature than that provided by the [O III] lines. We 
consider the first effect to be a minor one and the second to be a mayor one but 
both should be estimated for each object. 

From constant density chemicaly homogeneous models computed with 
CLOUDY we estimate that the maximum temperature that should be used to 
determine the helium abundance should be 5% smaller than T e (0 III). Moreover, 
if there is additional energy injected to the H II region T e (He II) should be even 
smaller. 

Luridiana, Peimbert, & Leitherer (1999) produced a detailed photoionized 
model of NGC 2363. For the slit used by Izotov, Thuan, & Lipovetsky (1997) 
they find an ICF(He) = 0.993; moreover they also find that the T e (0 III) 
predicted by the model is considerably smaller than observed. From the data 
of Izotov et al. (1997) for NGC 2363, adopting a T e (He II) 10% smaller than 
T e (0 III) and AY/AZ = 1.9 ± 0.5 we find that Y p = 0.234 ± 0.006. 

Similarly, Stasinska & Schaerer (1999) produced a detailed model of I Zw 18 
and find that the photoionized model predicts a T e (0 III) value 15% smaller than 
observed, on the other hand their model predicts an ICF(He) = 1.00. From the 
observations of A A 5876 and 6678 by Izotov et al. (1999) of I Zw 18, and adopting 
a T e (He II) 10% smaller than T e (0 III) we obtain Y p = 0.237 ± 0.007; for a 
T e (He II) 15% smaller than T e (0 III) we obtain Y v = 0.234 ±0.007, both results 
in good agreement with our determination based on the SMC. Further discussion 
of these issues is presented elsewhere (Peimbert, Peimbert, & Luridiana 2000b). 

The primordial helium abundance by mass of 0.2345 ± 0.0030(l<r) — based 
on the SMC — combined with the computations by Copi, Schramm, & Turner 
(1995) for three light neutrino species implies that, at the 95 percent confidence 
level, f^/i 2 is in the 0.0046 to 0.0103 range. For h = 0.65 the Y p value corresponds 
to 0.011 < fife < 0.024, a value considerably smaller than that derived from the 
pregalactic deuterium abundance, D p , determined by Buries & Tytler (1998) 



The Magellanic Clouds and the Primordial Helium Abundance 



9 



that corresponds to 0.040 < fi& < 0.050 for h = 0.65, but in very good agreement 
with the observational estimate of the global budget of baryons by Fukugita, 
Hogan, k Peebles (1998) who find 0.007 < Q b < 0.038 for h = 0.65. The 
discrepancy between Y p and D p needs to be studied further. 

To increase the accuracy of the Y p determinations we need observations of 
very high quality of as many He I lines as possible to derive T e (He II), iV e (He II), 
and t(3889) self-consistently. We also need observations with high spatial reso- 
lution to estimate the ICF(He) along different lines of sight. 

It is a pleasure to acknowledge several fruitful discussions on this subject 
with: L. Carigi, V. Luridiana, B. E. J. Pagel, M. T. Ruiz, E. Skillman, G. 
Steigman, S. Torres- Peimbert, and S. Viegas. 
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